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An emission tomography of laser-induced plasmas employed in the laser induced 
breakdown spectroscopy (LIBS) requires long signal integration times during which 
the plasma cannot be considered stationary. To reduce the integration time, it is 
proposed to measure a plasma absorption in parallel rays with an aperture that 
collects light coming from large fractions of the plasma plume at each aperture 
position. The needed spatial resolution is achieved by a special numerical data 
processing. Another advantage of the proposed procedure is that inexpensive linear 
CCD or non-discrete (PMT, photodiode) detectors can be used instead of costly 
2-dimensional detectors. 

1 Large apertures and the Abel inversion 

It is assumed that a plasma plume created by a laser ablation is axially symmetric (the 
symmetry axis coincides with the laser ray). Let the coordinate system be set so that 
the symmetry axis is the z— axis. Then the emissivity is a function e = e(r,z,t,u) 
where r = (y 2 + x 2 ) 1 / 2 , t is the time, and v is the frequency. A measured quantity is the 
intensity I(y, z, t, v) of light per unit time and unit frequency along the rays through an 
infinitesimal area element AA = AyAz centered at the point (0, y, z) that are in a narrow 
solid angle AVL about the line parallel to the 



It is well-known that Eq.([T]) can be solved for e(r, z, t) by the Abel inversion [HE]- A 
plasma plume has a finite size that determines a cut-off of the infinite integration limits in 
([1]). Equation ([1]) is valid only for infinitesimal A A and AQ and so is the approximation 
of the parallel rays that allows for the subsequent Abel inversion. In typical LIBS experi- 
ments, the smallness of AA is provided by a narrow spectrometer slit and small pixel size 
of a detector, while AQ is small due to a low acceptance angle of a spectrometer (high 
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/—numbers). To compensate for small factors AA and AQ, a long integration time T 
is used to collect enough energy E = Jq dtl(y, z, t, v) of photons emitted by the plasma 
along parallel rays. A typical integration time in LIBS plasma experiments is T ~ 1 fis 
or even higher. If the plasma is not stationary during the time T, the collected data 
correspond to a time-averaged intensity and the reconstructed emissivity may not be ac- 
curate. This is especially relevant for studying elemental contents of plasmas expanding 
into an ambient gas because of percolation processes at a rapidly moving plasma-gas in- 
terface. Numerical simulations of LIBS plasma dynamics show that this is often the case, 
especially for earlier stages of the plasma evolution (see, e.g., [U H]). 

Here it is proposed to increase an aperture through which parallel rays are collected 
in order to reduce the integration time so that the observed intensity is: 

ry+A poo ry+A 

I E (y) = AzAQ du dxe(Vu 2 + x 2 ) = / dul{u) . (2) 

Jy-A J-co ^ ' Jy—A 

where the dependence on z, t, and v is now suppressed as it is irrelevant for the discus- 
sion. Equation (j2J) states that all photons traveling 'parallel to the x— axis and coming 
through a rectangular aperture of an infinitesimal height Az and finite width 2A and 
centered at a distance y from the x— axis are collected. An experimental setup to which 
Eq. (T5]) applies is discussed below (see Fig. [TJ). In the limit A = Ay/2 — > 0, Eq. (2) turns 
into Eq.flTJ. So, the signal can be amplified roughly by the factor of 2 A/ Ay. A spatial 
resolution Ay in a typical experimental setup [5], El [7] is determined by a pixel size of a 
camera used to register I(y), Ay ~ 0.02 mm. If the aperture width is taken to be of a 
typical size of a LIBS plasma plume, 2A ~ 2 mm, up to two orders in magnitude of the 
signal strength can be gained. How is then the spatial resolution restored? It is proved 
below that, given a function Ie{v), the function I(y) can be uniquely determined from <^ 
for any A, provided I(y) satisfies the condition that I(y) vanishes for \y\ > D for some 
D, which is always the case for a physical I(y) because a plasma plume has a finite size. 
Thus, the proposed procedure entails: 

(i) Taking intensity measurements of Ie(u) at a set of positions y = y n ; 

(ii) Reconstructing the intensities I(y n ) from the data Ie^Uti)', 
(Hi) The Abel inversion for I(y n ). 

A key issue for the data collection (i) is to separate parallel and angled rays passing 
through the aperture. Note that for a non-infinitesimal A there should be rays coming 
through the aperture at a finite angle to its normal and, hence, AQ can no longer be 
considered infinitesimal, thus, invalidating the approximation of parallel rays in (121) . This 
is indeed true if the plasma emission is collected by the slit. The problem can be avoided 
if the collimated light absorbed by the plasma is measured. An experimental setup to 
achieve this goal is presented in Section 2 (see Fig. [[]). 

An explicit algorithm to carry out (ii) is given in Section 3. The spatial resolution is 
determined by the difference A s = y n — y n -\ in positions of the center of the aperture. 
Numerical simulations with synthetic data representative for LIBS plasmas are presented 
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Figure 1: An experimental setup for plasma absorption measurements with a large aperture. A 
CW light source is in focus of a lens to create a flux of parallel rays. A plasma plume with the 
symmetry axis normal to the figure plane is illuminated by the parallel rays. The parallel flux 
across a rectangular aperture is focused onto a detector (a spectrometer slit or CCD camera). 
The aperture plane is normal to the flux. It is positioned at a distance large enough to neglect 
contributions of the plasma plume emission. The aperture geometry is depicted in the inset at 
the bottom of the figure. The height Az is small enough to neglect variations of the plasma 
absorption along the z— direction. The aperture width 2A can be of the size of the plasma 
plume or larger. The intensity in parallel rays is measured for different values of y, the distance 
between the plasma symmetry axis and the optical axis of the system. Measurements at a fixed 
value of y have no spatial resolution by themselves if A is large. The spatial resolution is restored 
by a numerical processing of data collected at different values of y. 



in Section 4. There is a variety of numerical methods available to carry out the Abel 
inversion (see, e.g., [8], and more recent applications to LIBS plasmas (6j[7]). So the part 
(Hi) will not be discussed here. 



2 Plasma absorption experiments 



The proposed hypothetical experimental setup is shown in Fig. [TJ A point light source is 
placed in a focus of lens so that after the lens the light propagates parallel to its optical 
axis. The light goes through a plasma plume. Some of the photons will be absorbed 
by the plasma and then re-emitted. However the re-emission occurs in the whole solid 
angle of Air, i.e., the absorbed light no longer contributes to the parallel flux. Since the 
intensity of the plasma radiation is inversely proportional to the squared distance from 
the plasma center, by placing the aperture at a distance large enough from the plasma 
center the contribution of the plasma radiation through the aperture becomes negligible 
as compared to that of the parallel light. The plasma emissivity, e(r), is determined by the 
absorption coefficient n(r) and the black body radiation function B(r), e(r) = n(r)B(r), 
provided, of course, that the plasma is at a local thermodynamic equilibrium (see, e.g., [9]). 
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So, by measuring the difference between the fluxes in parallel rays piercing through the 
aperture area AA = 2AzA with and without the plasma plume inserted, the contribution 
of the plasma emissivity expressed as the attenuation of the parallel flux across AA is 
obtained. The plasma plume symmetry axis is positioned at a distance y from the optical 
axis of the right lens (the lens behind the slit). The center of the aperture also lies on 
this optical axis. The parallel rays that came through the aperture are focused on the 
spectrometer slit. In this setting, the spectrometer-detector combination is only used to 
achieve the spectral resolution because all the light that came through the aperture of 
width 2 A is focused to a single point of the spectrometer slit. Thus, there is no spatial 
resolution in measurements at a fixed value of y. An analogy to the data collection 
in the aforementioned emission experiments, e.g. P, [7], would correspond to the case 
when the right lens is removed in Fig. [I] so that the light is collected by all pixels of a 
detector evenly distributed along the spectrometer slit and both the spatial and spectral 
resolutions are achieved. In contrast, here all the light is focused onto a single pixel by 
the right lens in Fig. [TJ Thus, the signal amplification is achieved at the price of loosing 
the spatial resolution in measurements. However, as noted before, the spatial resolution 
is restored by a numerical data processing. It is noteworthy that 2D CCD detectors are 
no longer required to spatially resolve the plasma radiation. Inexpensive linear detectors 
(if the spectral dimension is to be retained) or photomultiplier tubes (PMT) can be used 
instead. 

The measurements of lE{y) are taken at different positions of the center of the aperture 
y = y n = nA s , n = 0, ±1, ±2, ±iV, where A s is the slit center displacement relative 
to the plasma plume symmetry axis (it is more convenient to create plasma plumes at 
different positions y n ). The step A s can be made as small as desired (or possible to 
achieve). The limiting positions y = ±NA S are chosen so that the studied plasma plume 
absorption does not contribute to the parallel flux through the aperture (i.e., NA S — A 
exceeds the plasma plume radius). 

The light source for absorption measurements must be bright enough, i.e., comparable 
to the plasma own spectral emission, otherwise the gain by the large aperture is rendered 
useless by a flux weaker than that of the plasma emission. Possible solutions of this 
problem are as follows. First, another laser-induced plasma plume can be used as the point 
light source. It is just as bright as the studied plasma plume, has a broadband spectrum, 
and a short pulse duration suitable for time-resolved measurements. Plasma plumes as the 
light source for absorption measurements have been used in experiments reported in [10] . 
Second, super- continuum white light lasers [11] with a broadband output spectrum can 
be used as a light source for absorption measurements. Finally, if the spectral resolution 
is not relevant, a single frequency laser can be used to study plasma emissivity at a 
particular frequency that coincides with one of the plasma emissivity spectral lines [T2] . 
The problem of an insufficient spectral brightness is also resolved in this case. 
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3 The reconstruction algorithm 

The integral equation (jSJ) is a linear equation. So its solution is the sum of a general 
solution Io(y) of the homogeneous equation (when I E = 0) and a particular solution of 
the non- homogeneous equation (j2J). It is assumed that I E {y) cannot be a non- vanishing 
constant (a plasma plume has a finite size). Differentiating Eq. (j2|) with respect to y and 
shifting the argument y — > y — A, one infers: 

I(y) = I(y-2A) + I' E (y-A). (3) 

If Ie = 0, then I(y) = Io(y) = Io{y— 2A) is a periodic function that is a linear combination 
of exp(iiiky/A) for k = ±1,±2, ... (the case k = is not possible for I E = in Eq. 
(T5])). By physical boundary conditions, I E {y) — I(y) — for all |?/| > D and some .D. 
No periodic Io(y) satisfies this condition. Thus, under this boundary condition, I(y) is 
uniquely recovered from I E {y). The solution can be written in the form: 

1 oo 

Hy) = ~ £{^(l/ " A ( 2 <? + !)) - r E(y + A ( 2 <? + W 

/ 9=0 

It can be verified by its substitution into the right side of Eq. (J2]). To understand the 
structure of the solution (jlj), it might be instructive to work out a simple analytic example 
of reconstructing I(y) = A = const for \y\ < A and I(y) = otherwise. The function 
I E (y) is easy to compute by Eq. (J2D, Ie(v) = A(2A - \y\) if \y\ < 2A and = 

otherwise. Then Ie(v) can be substituted into the right side of Eq. (TjJ to see how I(y) 
is recovered through cancellations in the sum. In this simple case, the function I(y) is 
reconstructed in the interval \y\ < A by the first term q = in the sum (jl]). The other 
terms vanish for \y\ < A and are needed to make I(y) = for \y\ > A. So, for a plasma 
plume of a radius R, only first i?/A terms are needed to reconstruct I(y) in the interval 

\y\<R- 

Put A = mA s for some integer m. Suppose that the data function /^(y) is taken at 
the grid points y n = nA s where n = —N, —N + 1, N — 1, N so that /^(y) = for all 
y < V-n and y > y^. Consequently, this implies that I(y) = for all y < y^^+m and 
y > yN-m- The number is chosen so that (N — m)A s exceeds the plasma plume radius. 
In practice, the data collection should simply start at \y\ large enough to see no signal 
from the plasma. Then y is changed with the step A s until the signal vanishes again. 
Equation ([3]) becomes the recurrence relation: 

With the boundary conditions imposed on the data set, it follows from the first 2m 
relations in (jSJ) that 

I{y n ) = I' E {y n -m) , n = -N + m+l,...,-N + 3m. (6) 
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Figure 2: Left panel: The emissivity as a function of r = \y\ (the horizontal axis) as described 
in the text (arbitrary units). The plasma radius is roughly R = 10. Middle panel: The intensity 
in parallel rays I(y) calculated by Eq. (JTJ) with the factor AAAfi omitted. It is normalized 
so that its maximal value is one. The (red) dots indicate the values I(y n ) reconstructed from 
lEiVn) for A = i? by the data processing algorithm of Section 3. Right panel: The intensities 
Ie(u) for three different aperture widths. The upper, middle, and lower curves are the graphs 
of I E (y) computed by Eq. ^ with A = R, A = R/2, and A = U/10, respectively. 

The values I(y n ) for n > —N + 3m are then determined recursively by ([5]) because the 
preceding values I{y n -2m) needed to initiate the recurrence relation ([5]) are now known 
from ()6]). The remaining problem is to calculate F E (y n ) fro m the set /e(|/ n ). 

Let lE{k) = J dylE{y)e lky be the Fourier transform of /e(|/)- The Fourier transform 
of F E (y) is iklE{k). Thus the values F E (y n ) can be found by a discrete Fourier transform: 

~- — -pi — 1 

(K) ^ I' E {y n ) , (7) 

where F implies taking the fast Fourier transform of the data set lE{yn) to obtain the 
set ijs(A; n ), and F~ l is the inverse of F. For analytic functions, the accuracy of a numer- 
ical differentiation by the spectral (Fourier) method is superior to any finite differencing 
method [13J. The spatial resolution of the reconstructed intensities I(y n ) is determined 
by A s . It could even be better than in conventional settings where spectrometers with 
a pixel detector are used for the spatial resolution, provided it is possible to achieve A s 
smaller than a spectrometer pixel size. 

4 Numerical simulations with synthetic data 

The emissivity is taken in the form e(r) = 6q exp(— r 2 /2a 2 )(a — cos(6r)) which resembles 
the emissivity profile typically observed in LIBS plasmas [7]. The parameters are chosen 
as Eq = 1, a = 3, a = 1.4, and b = 0.5 (arbitrary units). The graph of e(r) is shown 
in the left panel of Fig. [2j As seen from the figure, the plasma plume radius can be set 
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as R = 10 in these units. The intensity I(y) in parallel rays is calculated by means of 
Eq. (EE]) (the geometrical factor AAAQ is omitted). It is shown as the solid (blue) curve 
in the middle panel of Fig. |2j The function I E {y) calculated by means of Eq. (j2J) is 
shown in the right panel of Fig. (j2j). The top curve corresponds to the aperture width 
2A = 2R. For the middle curve, the aperture width is reduced twice A = R/2, and 
the bottom curve corresponds to A = R/10. In the conventional experimental settings, 
about 10 to 20 measurement across the plasma diameter are taken. So the bottom curve 
is representative for this case. A few remarks are in order. First, a significant gain in the 
signal amplitude is evident for the wide aperture. Second, the function Ie(v) has support 
twice as wide as the plasma diameter if 2 A = 2R because for R < \y\ < 2R no parallel ray 
comes through the plasma, but the aperture centered at such y can still capture parallel 
rays coming through the plasma. So the spatial resolution seems to be lost in the data 
Ie{u)- Third, the maximal amplitude of Ie{v) cannot increase any further if 2 A > 2R. 
So, the aperture width equal to the plasma size gives the maximal amplification of the 
signal. 

To illustrate the restoration of the spatial resolution, the function I E {y) is sampled 
at y = y n = nA s where A s = R/20. For the synthetic data I E (yn) when 2A = 2R, a 
discrete Fourier transform is carried out to compute I' E (y n ) as specified in (JjTJ). The Fast 
Fourier Transform Matlab package has been used for this purpose. The values I(y n ) are 
recovered by means of the recurrence relation (jHJ) with the initial conditions (JHJ). They 
are shown by the (red) dots in the middle panel of Fig. [2j The dots lie on the graph of 
I(y), i.e., the reconstruction is highly accurate. Naturally, the accuracy of the emissivity 
reconstruction from the data I(y n ) would then be fully determined by the accuracy of a 
particular Abel inversion algorithm [8]. 

The proposed scheme can also be used to improve the signal to noise ratio when the 
integration time is fixed. In a typical LIBS experiment, the intensity data are collected 
from repeatedly created plasma plumes. The power of a laser pulse used for ablation 
varies from pulse to pulse, a material surface on which the ablation process takes place 
is not ideal, etc. All these effects generate a noise in e(r) and, hence, in I(y). As an 
example, a white noise r)(y) has been added to I(y) — » I v (y) = I(y) + T}(y). The values of 
Ie{v) have been calculated by Eq. (TSJ) with I n (y) for 50 samples of rj(y), and the result has 
been averaged. Then the reconstruction algorithm has been applied to recover I(y) from 
the "noisy" data Ie(v)- Figure [3] shows the graphs of the reconstructed functions I(y) 
(dotted curves) for A = R/10 (left panel), A = R/2 (middle panel), and A = R (right 
panel). The solid (red) curve in each panel is the actual intensity I(y) corresponding to 
the emissivity shown in the left panel of Fig. [2J In contrast to I(y) depicted in the middle 
panel of Fig. &\ here I(y) is given without the normalization on its maximal value, i.e., as 
defined by Eq. (p]) without the factor AAAQ. The signal-to-noise ratio S = I m /cr, where 
a is the standard deviation of the white noise and I m is the averaged signal amplitude 
near its maximum, is S = 0.7, S = 2.9, and S = 3.1 for the above three values of A, 
respectively. The increase of the signal-to-noise ratio in the reconstruction of I(y) is 
self-evident. 
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Figure 3: Noise effects on the reconstruction algorithm. The solid (red) curves on each panel 
is the graph of I(y) as in the middle panel of Fig. [2] but without normalization on its maximal 
value. The dotted curves are the graphs of the reconstructed I(y) from the "noisy" data Ie{v) 
as explained in the text for A = i?/10 (left panel), A = R/2 (middle panel), and A = R (right 
panel). The signal-to-noise ratio is S = 0.7, S = 2.9, and S = 3.1, respectively, in these three 
cases. 
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